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Abstract 

In the next-to-minimal supersymmetric model (NMSSM) a light CP-odd Higgs boson is so far 
allowed by current experiments, which, together with a large tan /?, may greatly enhance the rare 
dileptonic decays B Xgl'^l" and Bg — > We examine these decays paying special attention 

to the new operator allowed by the light CP-odd Higgs boson. We find that in the parameter space 
allowed by current experiments like LEP II and b 57, the branching ratios of these rare decays 
can be greatly enhanced and thus the existing experimental data on B ^ Xsfx^n^ can further 
stringently constrain the parameter space (especially the region with a super-light CP-odd Higgs 
boson and large tan /?) . In the surviving parameter space we give the predictions for other dileptonic 
decay branching ratios and also show the results for the forward-backward asymmetry. 

PACS numbers: 14.80. Cp, 13.85. Qk,12.60.Jv 
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I. INTRODUCTION 



Recently some non-minimal supersymmetric models such as the next-to-minimal super- 
symmetric model (NMSSM) have attracted much attention [l| since these models can solve 
the /i-problem and alleviate the little hierarchy. In the NMSSM, for example, the /i-term 
in the superpotential is forbidden by imposing a discrete Z3 symmetry and instead it is 
generated through the coupling between the two Higgs doublets and a newly introduced 
gauge singlet scalar which develops a vacuum expectation value of the order of the SUSY 
breaking scale. In this way, the /i parameter at the weak scale can be naturally explained. 
The NMSSM can ameliorate the little hierarchy by either tuning the parameters to enhance 
the theoretical upper bound for the mass of the lightest CP-even Higgs boson or relaxing 
the LEP2 bound of 114 GeV through allowing for a light CP-odd Higgs boson (Ai) with 
mass below 2mb 

It is interesting to note that in the NMSSM the lightness of such a CP-odd Higgs boson 
can be naturally predicted in the enlarged parameter space, and is also allowed by the LEP 
II data [l|. This light Higgs boson can not only alleviate the littje hierarchy, but also can 
help to explain the observed anomaly in the decay S+ p/x+yU^ [3]. On the other hand, if 
a Higgs boson is indeed so light, its effects in some low energy processes may be sizable and 
thus are necessary to eheek fl B. For a super hght A, the deeay , ^ A,., .s opeu and an 
analysis has been performed in [41] (note that as analysed in [41], a small CP-odd Higgs mass 
is only protected from RGE-effects in the limit of large tan/?). In this work we consider the 
full possible mass range of (heavy, intermediately heavy and light) and check the NMSSM 
effects in the rare B meson dileptonic decays B — > Xs£~^£~ and Bg —>■ £"'"£"7 jo]. 

These rare dileptonic decays are induced by the flavor-changing neutral-current (FCNC) 
b ^ s transition and are of special interest due to their relative cleanness and high sensitivity 
to new physics. In the Standard Model (SM) such FCNC processes are suppressed and have 
very small branching ratios [7] but can be greatly enhanced in some new physics models 

iil,Q,Q- 

Since experimental data on B Xgfi'^fi is available and the future 



0,3,0, 



LHCb or super B factory will further scrutinize B meson decays, these dileptonic decays 
serve as a good probe of new physics. 

In supersymmetric models these dileptonic decays can be drastically enhanced by large 
tan/3 since both the 6 — > s transition loops (such as the charged Higgsino loops) and the 
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Higgs couplings in the Higgs-propagated diagrams are proportional to tan/5. It has been 
shown (see e.g. 12j) that in the minimal supersymmetric model (MSSM) great enhance- 
ments are possible for these decays. In the context of NMSSM, in addition to the tan/5 
enhancement, the presence of a light CP-odd Higgs boson could further enhance these dilep- 
tonic decays. Thus, the parameter space, especially the region with a super light CP-odd 
Higgs boson and a very large tan/?, is constrained by the existing data on B XsyU+yU". 
In our analysis we will examine the NMSSM effects in these dileptonic decays by scanning 
over the parameter space allowed by the LEPII experiments and the data on 6 — 57. We 
will show the 2a constraints from B Xsfi^fi^ on the parameter space and then give the 
predictions for other dileptonic decay branching ratios and forward-backward asymmetry. 

A key point in our calculations is the presence of a new operator due to the light CP-odd 
Higgs boson. In contrast to the MSSM, where all Higgs bosons and sparticles are so heavy 
that they can be integrated out at the weak scale, the light CP-odd Higgs boson Ai in the 
NMSSM cannot be integrated out at the weak scale and thus a new operator Oa describing 
the interaction Aibs must be treated carefully. 

This work is organized as follows. In Sec. [Tll a brief description of the NMSSM is 
presented. In Sec. IIIII we calculate the Wilson coefficients paying special attention to the 
new operator Oa- In Sec. |IV] we evaluate the NMSSM effects on the dileptonic decay 
branching ratios and the forward-backward (FB) asymmetry, and present some numerical 
results. The conclusion is given in Sec. |V] and analytical expressions from our calculations 
are presented in the Appendix. 

II. A BRIEF DESCRIPTION OF NMSSM 

In the NMSSM a singlet Higgs superfield S is introduced. A discrete Z3 symmetry is 
imposed and thus only the cubic and trilinear terms are allowed in the superpotential. The 
Higgs terms in the superpotential are then given by 

XSHu-Hd+'^S^. (1) 

Note that there is no explicit /i-term and an effective /^-parameter is generated when the 
scalar component (S) of S develops a vacuum expectation value s/\/2: 

^^eff = x{s) = ^. (2) 



The corresponding soft SUSY-breaking terms are given by 
The scalar Higgs potential is then given by 



D 



92 



\Hu?) 



(3) 

(4) 
(5) 



K 



y,oft = m',\HaV + mi\H^\' + mi\S\'+ [Ay^XSH^ ■ H, + -A^S' + h.c.j , (6) 

where = gf + with gi and g2 being the coupling constants of f/y(l) and SUl{2), 
respectively. 

The scalar fields are expanded as follows: 



^ [Vd + (Pd + i'^d) 
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The mass eigenstates can be obtained by unitary rotations 



^S = ^{s + a + tO . (7) 
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where -f^i,2,3 and Ai^2 are respectively the CP-even and CP-odd neutral Higgs bosons, G° 
and G~^ are Goldstone bosons, and H'^ is the charged Higgs boson. It is clear that the 
charged Higgs sector is the same as in the MSSM, while the neutral Higgs sector contains 
one more CP-even and one more CP-odd Higgs boson. U"^ and are given by 
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tan /3 s tan fS ^ J 



where Gx = cosX and Sx = sinX {X = 9a,6h)- The mixing angles are given by [l^ 

1 \ . 2\'^vs 



2 s tan P 
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with V ~ 246 GeV and 
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V2Ax =F 2ks 
V2Ax + Ks 
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The Lagrangian of the Higgs couphngs to quarks for a large tan (3 are given by 

^A,M = -^1^ tan/3A) J75d, (12) 

L^Aiuu = -It, ; — ^ — — M u-f5U, (13) 

>C(^^„if„if3)<i<i = -f^ ((<^^ - (^e + l5+Ce)H2Mn(3H^) dd, (14) 

^(^,if.if3).« = t^t^. (15) 

Here one can see that the couphng of the neutral CP-odd Higgs Ai with up-type quarks are 
suppressed by a large tan/3 and thus can be neglected in the large tan/3 limit. 

Since one more Higgs superfield 5" is introduced in the NMSSM, we have a new neutral 
Higgsino ■05- So the neutralino sector is composed of C/y(l) gaugino A^, SUl{2) gaugino A^, 
and the Higgsinos -0)^^, y/l^^ and ■05. The corresponding mass terms are given by 

= i\v{g2\^ - g^\^){cos - sin/JV'IJ - ^M^A^A^ - ^M^AU^ - -La^V^^.^V'^ 

(cos/^V^I^^s + s\n(3il)]jj)s) + -^'^^'^l + li-c- 
= -^(^°)^K,o^° + h.c., (16) 

where 

(V'°)^ = (-^A^-^A^V'l,„V'k,V'5). (17) 

The neutralinos are obtained by the unitary rotation ■^^ = (Ziv)ijX°, where Zjq diagonalizes 
the mass matrix Y^o. 

Similar to the charged Higgs sector, the chargino sector of the NMSSM is the same as in 
the MSSM with ji replaced by //e//- The chargino masses are obtained by the diagonahzation 
of the mass matrix with two unitary matrices Z_ and Z^: 

rj. I —Mo -\/2mwsin/? \ 
M^c = {Z. f I ^ ^ ^ U+. (18) 

yV^mw cos P -Mf ) 

III. CALCULATIONS OF WILSON COEFFICIENTS 

In our calculations we consider the flavor mixing between h and s, which make contribu- 
tions to the dileptonic B meson decays through gluino or neutrahno loops. Following the 



analysis in 15| , we assume the flavors are diagonal at tree level and the mixings are induced 



at loop level. Such mixings can be parameterized by a small mixing parameter ei which is 
dependent on some soft-breaking mass parameters [l^. In our numerical calculations we 
input ei = 0.1 for illustration. We perform the calculations in the Feynman gauge and thus 
the Goldstone bosons will be involved in the loop diagrams. 



3iggs boson Ai, 



l2| ( using the 



Since in the NMSSM the lighter mass eigenstate, the CP-odd neutral 
can be rather light, with a mass ranging from 100 MeV to the weak scale 
package NMHDECAY [ij], we checked that such a light CP-odd Higgs boson Ai is indeed 
allowed by the LEP II data), in our calculations we pay special attention to this wide mass 
range of Ai and discriminate three cases. 

(i) Case A: Heavy Ai. 

For a heavy Ai, around weak scale, we integrate it out together with the other heavy 
particles (Higgs bosons, top quark, and Z bosons, sparticles) at weak scale to ob- 
tain the Wilson coefficients. The effective Hamiltonian describing b sl^l^ transition 
reads 

10 2 
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(19) 



.1=1 i=l 

where Oi and Qi are operators listed in [sl, Q, and Ci and Cq. are respectively 
their Wilson coefficients, and fir is the renormalization scale. Note that the most 
general Hamiltonian in low-energy supersymmetry also contains the operators O'^ and 
Q'j which respectively are the flipped chirality partners of Oi and Qi. However, they 
give negligible contributions and thus are not considered in the flnal discussion of 



physical quantities 



17| 



In this case, there are no new operators and we only need to calculate the NMSSM 
contributions to these coefficients Ci and Cq. at the scale of ^m^r. For the processes in 
our analysis, only Cy^g^io and Cq^ j are relevant. The NMSSM contributions to 6*7,9,10 
are the same as in the MSSM, which are computed in [l^. For Cg^a NMSSM 
contributions are different from the MSSM contributions |l^ and thus need to be 
calculated here. The Feynman diagrams we need to calculate are shown in Fig.l, 
where the loops respectively involve the charged Higgs bosons, charginos, gluinos and 
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neutralinos. From the calculations of these diagrams we obtain the Wilson coefficients 
at mw scale, which are presented in the Appendix. 
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FIG. 1: The Feynman diagrams which give the dominate contributions to Cq.^^- (^"C) charged 
Higgs loops, (d-g) chargino loops, (h-k) gluino and neutralino loops. 

For the calculation of the dileptonic B meson decays we need to know the Wilson 
coefficients at the rub scale, which can be obtained from the running of the coefficients 



at the mw scale down to the sea 
dimension which can be found in 



12 . 



e. Such a running is governed by the anomalous 



(ii) Case B: Ai with an intermediate mass m^, ^ m^^ <^ 

At the m-w scale we retain Ai as an active field and thus we have a new operator 
Oa- After integrating out all heavy particles at the mw scale we obtain the Wilson 
coefficients including Ca for Oa- Then we work out the anomalous dimensions and 
run the Wilson coefficients from the mw scale to the m^^ scale. At the m^^ scale we 
integrate out Ai, which makes an additional contribution to CQ.(m^J. Finally, we 
run all the Wilson coefficients from the m^i scale to the mb scale. 

The new operator Oa at the mw scale takes the form 
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A = I 



;mbmwslb%Ai. 



(20) 
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From the calculations of the corresponding diagrams in Fig. 1 we obtain the Wilson co- 
efficient CA{fnw)- It is composed of the charged Higgs loop contribution, the chargino 
loop contribution, the neutralino loop contribution and the gluino loop contribution, 
whose analytic expressions are given in the Appendix. 

For the running of the Wilson coefficients including Ca from the m^y scale to the 
rriAi scale we work out the anomalous dimensions by calculating one-loop diagrams 
with operator insertions. We find that all the Wilson coefficients in Ea. ffTOl) run in 

n 

the same way as in the MSSM [12[, while the new coefficient Ca is not changed, i.e., 

CAiniAr) = CA{mw)- 

When integrating out Ai at the m^^ scale we find it gives a contribution ACQ^{mA) 
to the operator Q2 

ACQMA.) = -^-l'^CA{mA^ (21) 

Finally, for the running of the Wilson coefficients from the scale to the scale 
the anomalous dimensions are the same as in the MSSM [l^ . 

(iii) Case C: Super light Ai with mass ttiai < mf,. 

In this case we retain Ai as an active field in the entire analysis. At the rriw scale 
we integrate out all heavy particles and obtain the Wilson coefficients including Ca- 
Then we run the coefficients down to the ruf, scale. At the rrif, scale the effects of Oa 
are represented by a change in Cq^ 

ACq, mfe = — — —CA{mb), (22) 

2 sp^ -m\^ +imAj^A^ 

where p is the momentum transfer and F^^ is the total width of Ai. Since Ai can be 

on-shell in this case, the effects of Ai can be sizable even if without tan (3 enhancement. 

Note that the chargino loop contributions to Ca were also calculated in {4] where the corre- 
sponding diagrams induced by the Ai-squark-squark vertex are neglected since the author 
considered the large tan/3 limit. In our numerical calculations we used the full results by 
keeping all terms and thus we also included the diagrams induced by the Ai-squark-squark 
coupling although they contain no leading tan (3 terms. Except for the case of a large tan /3, 
such diagrams induced by the Ai-squark-squark coupling should be included since the Ai- 
squark-squark coupling can arise from the F-term of the superpotential and not suppressed 
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by the singleness of Ai. We checked that in the large tan (3 limit we can reproduce the 
analytical result given in ^\ for the chargino-loop contributions. 



IV. DILEPTONIC B-MESON DECAYS IN NMSSM 

With the effective Hamiltonian and the running of the Wilson coefficients presented in the 
preceding section we calculate the inclusive decays B Xs£~^i~ and their forward-backward 
(FB) asymmetry, as well as the exclusive decays Bs i^i~'~f. The formulas in terms of the 



Wilson coefficients can be found in 



121. 



Note that our supersymmetric contributions to the Wilson coefficients are given at one- 
loop level (next-to-leading order), while the SM contributions are known at two- loop level 
(next-to-next-to- leading order) [t!]. In our numerical calculations we consider the one- loop 
results for the NMSSM, while for the SM we also include the two-loop results. 

For the inclusive decays B Xsi^i~ we exclude the resonances J/\I' and contributions 
by using the same cuts as in the experiments 13], i.e., the invariant dilepton mass in the 
ranges 

(2m,, 2.75 GeV) © (3.3 GeV, 3.39 GeV) © (3.84 GeV, mt), (23) 

so that our results can be compared with the experimental measurements. 

For exclusive decays Bg 7 we follow jiol . and consider the photon in Bg — »• 

i~^i~'~f as a hard photon by imposing a cut on the photon energy E^, which means that the 
radiated photon can be detected in the experiments. This cut requires > 6 mBj'2 with 
6 = 0.02. (Note that for a soft photon both processes Bg — > and Bg must 



be considered together and in this case the infrared singu 



ar terms in Bg — > i^i 7 are be 
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) 



cancelled by the 0{aem) virtual corrections in Bg £~^i 

In our numerical calculations we perform a scan over the NMSSM parameter space 

2 < tan /? < 30, - 500 GeV < fx^fj < 500 GeV, 
-1 < A < 1, - 1 < K < 1, 

-50 GeV <Ax<50 GeV, - 10 GeV < < 10 GeV (24) 

with fixed parameters for the sfermion and gaugino sector (500 GeV for all sfermions and 
the gluino, and 200 GeV and 100 GeV for SU{2) and U{1) gaugino masses M2 and Mi, 
respectively). In our scan we consider the following constraints: 



(1) The LEP2 constraints by using the package NMHDECAY [16!]. 

(2) The constraints from B —>■ Xsj which stringently constrain the effective coefficient 
Cj^-^ . For the experimental result we use the world average value 20] 



Br{B ^ X,7)|e.p = (3.55 ± 0.24^^;^^ ± 0.03) x 10 



(25) 



(3) The constraints from Bs 



which constrain the Wilson coefficient Ca of the 



light pseudoscalar operator [j, ll2|. The experimental result is given by 21] 

BriBs AtV") < 1.5 X 10"^ (90% C.L.). (26) 

Among the relevant dileptonic decays experiment data is only available for Br{B — *■ 
Xsfi~^fi~), which is given by 21] 



Br{B Xsfi^fi-) = (4.3 ± 1.2) x 10" 



(27) 



In displaying our numerical results we will show this bound and use it to constrain the 
parameter space. For other dileptonic decay branching ratios, with no experiment data 
available, we will compare the NMSSM predictions with the SM values given by 



Br{B X,r+r") = 4.43 x 10" 



Br\B, 7r+r-) 



1.33 X 10" 
1.35 X 10^ 



(28) 
(29) 
(30) 



Note that the SM prediction for Br{B —>■ XgT^T^) was also given in [2^. But our result is 
different from theirs because it is very sensitive to the cuts around the resonances J/ ip and 
ip'. While our cuts are chosen as in Eg. (1231), w e cannot find the corresponding cuts used in 



22l |. We can easily reproduce the result in 



22l | by varying the cuts. 



In Figj2]we show the scatter plots of the branching ratio for B XsyU+yU^ versus tan/3. 
Here we present the results for the three cases: a super light Ai (m^^ <5GeV), an interme- 
diately heavy Ai (5 GeV< m^^ < 40 GeV) and a heavy Ai at rriw scale. In order to see 
how stringent the b s'~f constraints are we display the scatter plots with and without the 
b ^ s'y constraints. From this figure we make the following observations: (1) The branching 
ratio can be greatly enhanced by large tan/3. (2) b s'y constraints are quite stringent 
and can exclude a large part of the parameter space, typically with large tan /5. (3) In the 
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FIG. 2: Scatter plots for the branching ratio of 5 ^ Xsfj,~^^~ versus tan/3: the left panel is for a 
super light Ai {uia^ < 5GeV), the middle panel is for an intermediately heavy Ai (5 GeV < niA^ < 
40 GeV) and the right panel is for a heavy Ai at my\r scale. The dark (red) points are allowed by 
b — > S7, while the light (sky-blue) points are excluded by 6 — > 57. 

parameter space allowed by 6 — » S7 the decay B Xs^^ can still be greatly enhanced, 
especially in case of a super light Ai. The 2a experimental bound on i? ^ X^/i+yU^ can 
further exclude a large part of the parameter space. Almost no points with tan > 15 in 
the parameter space survive all the constraints. From the left panel of Fig|2] we see that 
some part of the parameter space with a super light Ai is still allowed by 6 — S7 and 
B ^ 

Let us take a look on the constraints from the process Bs — /^"'^/x" , whose branching ratio 
is given by [23 1 



Br{B, 



1.2 X 10"^ 

X C^n, + 
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+ 2 — —CiQ 



(31) 



We see that the contributions are from Ciq and Cg^ ^ ■ While the contribution from Cio is 
suppressed by the factor m^/niBs, ^^e contributions from Cq-^^ can be enhanced by large 
tan/3 (Cq^2 contain terms which are proportional to tan^/5, as shown in the Appendix). 
This feature can be seen from Figj3] in which we set aside the 6 — 57 constraints and 
illustrated the Bg — > fi'^fj^~ constraints. We see that, similar to the b ^ s'-f constraints, the 
Bs IJ'^/J'^ constraints are stringent for a large tan (3. If we impose the b constraints 
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which exclude a very large tan /3, then the further constraints from Bg — > are stringent 

only for the parameter space with a very light Ai. 
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FIG. 3: Same as FigJJl but with the Bs f^^fJ. constraints. The dark (red) points are allowed 
by Bs /^^/i , while the light (sky-blue) points are excluded by Bs fJ-^fJ-' ■ 



The results for other dileptonic decay branching ratios, for which no experiment data are 
available, are presented in FigHl To see how stringent the constraints from B Xs^^ ^~ 
are we display the scatter plots with and without such constraints (all the points satisfy 
b —>■ s'-f and Bs fi^fi^). 

From FigJUwe see that under the constraint from B Xsfi^fi^, the branching ratio of 
B —>■ XsT^T' does not deviate significantly from the SM value. The reason is that these 
two decays are highly correlated except for the contributions of Cq^ ^ and Ca which are 
dependent on the lepton mass. If the contributions of Cq^ ^ and Ca are dominant, then 
Br{B XsT^T^) should not be severely constrained by -B ^ XsfJ^^fi^ . As discussed below 
Eg. (1311) . the contributions of Cg^ ^ are important only for very large tan/5 which is 

not allowed by 6 ^ s'-f. As a result, the contributions of Cg^ j ^^"^ are not dominant and 
thus B XgT^T^ is highly correlated to S — Xsfi^fi~ . 

Note that for B X^r+r" it may be rather challenging to disentangle the NMSSM 
effects from the SM value in future experiments. One reason is that, as discussed above, the 
NMSSM effects are no longer so sizable under the constraint of -B — > Xsfi'^fi~. The other 
reason is that the SM prediction has its own uncertainty. If we consider the uncertainty of 
the input SM parameters, we can obtain the uncertainty (about 20% as found in [22]) of the 
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FIG. 4: Same as Fig|2l but for B XsT+t', Bs and Bs -fT+T'. The dark 

(red) points are allowed by -B — > X^/i^/i^, while the light (sky-blue) points are excluded by 
B^Xs^i+^i-. 



SM prediction. But in FigJH we did not show such an uncertainty of the SM value because 
for all the results, both NMSSM and SM, we used a same set of the SM parameters without 
allowing them to vary in the uncertainty range. Since the SM parameters are involved 
in both the NMSSM and SM values, all the results are subject to some uncertainty if we 
consider the uncertainty of the SM parameters. Of course, such uncertainties will deteriorate 
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the observability of the NMSSM effects. 



m^=4.212 




Q 0.2 0.4 0.6 0.8 

s 



FIG. 5: The forward-backward asymmetry oi B ^ Xsfi^fJ. versus s = /m^ {p^ is the invariant 
mass of and ^~ ). 

Finally, in Fig|5] we show the results for the forward-backward asymmetry m. B ^ 
Xsfi~^fi~ under the constraints from B Xsfi^fi~ . 

V. SUMMARY 

In the framework of the NMSSM we examined the rare dileptonic decays B Xsi^i~ 
and Bs i^i~'~f paying particular attention to the light CP-odd Higgs boson. We found 
that in the parameter space allowed by current experiments, such as LEP II and b — > 57, 
the branching ratios of these rare decays can be greatly enhanced and thus the experimental 
data on B Xg/i+yU^ further stringently constrains the parameter space, especially with 
a super light CP-odd Higgs boson and large tan/?. In the surviving parameter space we 
gave the NMSSM predictions for other unmeasured dileptonic decays which may hopefully 
be measured at the future LHCb or super B factory. 
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APPENDIX A: WILSON COEFFICIENTS 

The Wilson coefficients Cj, Cg and Cio in the NMSSM are the same as in the MSSM 
[isl . Here we present the new coefficient Ca, and Cq^ and Cq^, whose predictions in the 
NMSSM are different from the MSSM. We checked that we can analytically reproduce the 
"HQ (However, the NMSSM results can not explicitly reduce to the 



MSSM results 



12 



MSSM results by simply dropping out the singlet S (say setting X = k = 0) because the 
/i-term is generated by S). 

Although in our numerical calculations we used the complete results by keeping all terms, 
here, for simplicity, we only present the terms which can be enhanced by large tan f3. At the 
mw scale each Wilson coefficient is composed of the charged Higgs loop contribution from 
Fig.l(a-c), the chargino loop contribution from Fig.l(d-g), the neutralino and gluino loop 
contribution from Fig.l(h-k). 

For the charged Higgs contributions: 

Cf* = -^^tan/3Fi{xH±t,xwt) (Al) 
g2mw 



/id 



-.2 



-UalUalFl{xH±t,Xwt) + ^ ^ Fi{XtH± , Xtw] 



(A2) 



2 



2 2 

H 2 2^-^! (^tH± , Xtw) 



(A3) 



Note that although H^W^Ai vertex has 1/tan/? suppression by the singleness of Ai, 
H^G^Ai vertex which comes from the soft term A\\SH}^H^ does not have such 1/ tan/5 
suppression. Thus the contribution to Ca from the loop involving and with H^G^Ai 
coupling (in the Feynman gauge) is proportional to tan f3. 
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For the chargino contributions: 



(^tVixP ' ^xf xp) - ^nj^-^ (^tVixF ' ^xf xp) I ' 

3 2 ' 



4miT mw cos p ■' ''-i^j 



-SijPi (^tVi4-i ' ^xf 4-1 ) 

tk-i 



<3:ij^2(a;tVixp'^x?xp) + ;;r^^«i'^i(^tVixp'^xfxp) 



m r 



^- 4'^4(^tVixf ' ^xPxf ' ^i>X,^ 

- + ^3(^f,_iX?'^xPxf'^i>xP> 



r "i^c 



Oik 

92 



V2U^^Tf''mtm^c 

;5 —5liF-i (Xf f , X-Cf ) 

tk-i 



2 

H 2^<^ife 

^3 



^- ^+ ^4 ' ^xPxP ' ^^XP ^ 



~^c^- ^+ ^3(a;tV,;ef'^xPx?'^i>xP 



Note that here Cq^ .-, contain terms which can be 

V? fi in large tan /? hmit) . 



(A4) 



(A5) 



(A6) 

enhanced by tan^ j3 (the overaU factor 



J LlidL liClC l^Qj 2 L-LlliUCllli UClliiO V 

tan^ j3 multiphed by 1/ cos (3 gives tan 
For neutrahno contributions: 

rx° _ z tan/3 , / a^g* [" 1/2 ^ z x 

)l/2 



-^51'^2K_,^0,X^O^O)J+T3-T2(Z,J,/C')(X^, 

^ f \ 5-cot(5v 1/2. 



+T3i^'='r2(z,j,/c')(x^o^^,_^x^^^,_^ 

^^.V2,^(...'m(.,_.,o)} (A7) 
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' K,,K*4m^ cose ^2^2. 2^ M ^^^^ d,k'r2{z,j, k)P^{x 



92 



"^x? . 1 



al 2 



laik' y^j* rj-dl* rpk'2\ 



V2 



-iV_i'^x?V-i^ 



Xfc 



Co 



ti /la l,fe' = l J,fe=lm=l ^ 

Xj Xj 



al'^al"''X 



97-rA 73*;* 



m. 



Xfc 



!^r6(z, Zi-Z^^* -V,{i,k,j)\ F3(x^^_^,o , x,o,o , x,-^, 

Xfc ' 



For gluino contributions: 



CI 



tan pT^j^ 



1 nibmi IQgirng ^ 



+ 



-Fi(xr r ,X-r ) 



+ 



-Fi(x; 



19' 



bi-ibj-i'^fbj-i' 



bj-i 
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In the above expressions the constants and functions are defined by 

N'j = ^Z]j sin 9yj - cos 9yj, N'- = -Zjj sin 9yj + cos 6, 

Qai, = ^{u^,z^-^zl^ + u^,z}^zl^) - ^u^,z^Jzl^ 



= + ztjK) + V2Xiz'^z%' + z^z^ 

^^^^^-tB^^^^^K + ^N^l) - ^Kz'^zf + z'^zf)]} 

COo (7' 



w 



-V2KU^,iz'JZ%' + Z'JZ%') 

Q% = Q^k = \{u^2\-^{zj^N'; + Z^jK) + V2XiZ^Zf^ + Z'^Z%^)] 
+^ai[-^iZ%'N. + Z^X') + V2Xiz'jZj,' + Z'^Z%')]} 

-V2kU^,{z'Jz%' + z'JzI,') 

T^{^,k,J,l) = {T^n^- - 6,,8,,)Zf Z'f - . TifT^^'zfz'f 

yzmw sm p 

r,ii,j,k') = V2tan0^4^>i^*T^'3 + T^'^ 

V^JTivK COS p 

3m6 COS 6't„ 

r7(^,i,A;) = -^z^^T^^T^izl-'* zl^{zf sine^ + zf cose^)zf'' 

Zi COS ^ix) 

-2 sin 0^4^+^)™' Z^'^*] 

Tiaife A rpi3* rpk2 A*rpi2'' rpk3 rpaik a rpk2* rpii A*rpk3*rpi2 

As = ^{vdU^3 + sU^,) - AuU^2, A, = -^{vdU^s + sU^) + AuU^^ 
^5 = (1 - i sin' ^"')^a1, ^6 = \mw tan^ 

Z cos Uyj o o 

Ay = mwedi&ri' 0yjUa2, ^8 = + 2ed sin 6'^)[/J 

^ COS 

xlnx , 1 /xlnx y\ny\ 

F2{x,y) - , , 
X — y \ X — 1 y — 1 J 

/ N xlnx 18 , , , 

F3(x, y, z) = — + (X ^ y) + (x ^ z), 

(x — l)(x — y)(x — z) 



X — 1 X — y\x — 1 y — 1/ 

1 /x^ Inx lny\ 
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